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Abstract 



Problems posed by semirelativistic Hamiltonians of the form H = ^/rn?+p^ + V {r) 
are studied. It is shown that energy upper bounds can be constructed in terms of 
certain related Schrodinger operators; these bounds include free parameters which 
can be chosen optimally. 
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1 Introduction 



We study semirelativistic Hamiltonians H composed of the relativistically cor- 
rect expression Kip^) = \/ + , p = |p|, for the energy of a free particle 
of mass m and momentum p, and of a coordinate-dependent static interac- 
tion potential V{t\ r = |r|, which may be chosen arbitrarily, apart from the 
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constraint imposed on H that it be bounded from below: 



H = ^Jw?+p'^ + V{r). (1.1) 

The eigenvalue equation generated by this kind of Hamiltonian is usually 
called the spinless Salpeter equation. It arises as a well-defined approximation 
to the Bethe-Salpeter formalism for the description of bound states within 
(relativistic) quantum field theory [1] when it is assumed that the bound- 
state constituents interact instantaneously and propagate like free particles 
[2]. At the same time, H may be regarded as the simplest and perhaps most 
straightforward generalization of a (nonrelativistic) Schrodinger operator to- 
wards the incorporation of relativistic kinematics. For many potentials, this 
Hamiltonian can be shown [3] to be bounded below and essentially self-adjoint, 
and its spectrum can be defined variationally. For definiteness, we consider the 
corresponding eigenvalue problem in three spatial dimensions. 

In Sec. 2 we review the well-known tangential Schrodinger upper bounds which 
may be found [4] -[6] either by use of optimized operator inequalities or by the 
exploitation of the concavity of the Salpeter kinetic-energy X as a function 
of p^. The new Schrodinger bounds which are the principal concern of this 
paper are derived by considering operator difi'erences. We shall now illustrate 
the main idea by considering a nonrelativistic example. Suppose we wish to 
estimate the bottom of the spectrum of 

H = p^ + ar^- (3r^ ^ Hi - H2, (1.2) 

where a and /3 are positive and 

Hi = {l + uj)p^ + ar^ and H2 = ujp^ + /3r^ a; > 0. (1.3) 

Since Hi = H + H2, we conclude from the theorem of Weyl [7]- [9] that Ei > 
E + E2. We note in passing that for the ground-state energies discussed here, 
the Weyl inequality follows immediately by applying the exact normalized 
wave function ■^i corresponding to i/i as a trial function for the terms of the 
sum. Thus we have 

El = (V^i, Hi^i) = (V^i, H^i) + (V^i, H2^i) >E + E2. (1.4) 

Of course, we assume that the operator domains allow this. It remains to 
optimize the expression for E with respect to a; > 0. Thus we find in the 
example 

E < min [£'1(0;) — E2(oj)] — min 

where in three dimensions we have for the respective spectral bottoms 

+ / ^ = 64 f« 3.799673, and + ^ E ^ 62 ^ 3. (1.6) 
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The coupling dependence is found by the general scaling law 

p^ + v sgn(g) E{v) = E{1) ^;2/(2+9) (17) 

for pure powers ^ q > —2. For our problem and the special case a — (3 — 1, 
we find the result 

2.8345362 ^ E < ^ 2.85525 for cu = 0.818584. (1.8) 

By this reasoning we determine the energy of if = + — with error 
less than 0.74%. In Sec. 3 we shall show how this idea can be applied to the 
Salpeter eigenvalue problem. 



2 Tangential Schrodinger upper bounds 



The kinetic-energy term Kip^) = \/ + in the Hamiltonian H = K -^-V 
is a concave function of . Thus tangents to K generate upper Schrodinger 
operators HS^^ of the form 

ii-<7iW =a(t)p2 + 6(i) + l^(r), (2.1) 

where i > is the point of contact between the tangent ap^+b and K{p^). Ele- 
mentary analysis allows us to obtain the following formulas for the coefficients 
a{t) and h{t): 

zyrri'^ + t lym^ + t 

If an eigenvalue of the Schrodinger operator ap^ + V{r) is given by ^(a), then 
we have by the variational characterization of the discrete spectrum of H that 
the corresponding eigenvalue E oi H is bounded by the inequality 

E<rmn[£{a{t)) + b{t)]. (2.3) 

The minimum in this expression simply picks out the energy of the best upper 
tangential operator. We have shown earher [4] that these 'envelope bounds' 
are identical to those obtained by optimizing over the parameter /i in the 
upper bound for K{t) implied by the inequality ||A'(t) — > 0, namely 
K < {K'^ + /x^)/(2/i); the link between the two expressions for the bound is 
the parameter relation n = y/mF+t. 

The advantage of the tangential bound is its generality: it applies to each 
discrete eigenvalue that exists for the upper operator. For later comparison 
we consider three examples. We restrict our attention to the lowest eigenvalue, 
which is the subject of the difference upper bound discussed in Sec. 3 below. 
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2.1 The ultmrelativistic harmonic oscillator H — p + r^ {m — 0) 

By the spectral equivalence H = H = + r, we see that the exact energy 
is given by —zq, where zo is the first zero of Airy's function Ai(2;). That is 
to say E fa 2.3381074. In order to compute the envelope bound we may re- 
parametrize (2.2) in terms of s = a = 1/(2 Vt) and find b = l/(4s). The 
tangential Hamiltonian then becomes Ti = sp^ + l/(4s) + r^. The lowest 
eigenvalue of this operator is then given by £ {s) = 3 1 / (4s) . If we minimize 
£{s) with respect to s we find E < (9/2)6-^ ^ 2.47644. This is about 5.9% 
high. 



2.2 The semirelativistic harmonic oscillator H — -\-p^ -\- (m = 1) 

In this case the Hamiltonian H is equivalent to the Schrodinger operator 
H — p^ + vT+~r^ whose exact energy E is straightforward to find numerically 
and is given by E Ki 2.6640196. Meanwhile the tangential operator is given by 
TC = ap^ + b + r"^ and has lowest energy £{f) = 3^a{t) + b(t), where a{t) and 
b{t) are given by the formulas (2.2). A minimization of S{t) with respect to t 
yields the best upper bound E < S{3) = 11/4. This bound is about 3.2% high. 
As m is increased in the operator H = \/m^ + + r^, the problem spectrally 
(and monotonically) approaches the Schrodinger hmit B. — m-\-p^ /{2m) -|-r^, 
for which the envelope approximation is exact. 



2.3 The ultrarelativistic linear potential H — p + r {m — 0) 

This very symmetrical operator is truely non-local but yields easily to a varia- 
tional treatment in a Hermite basis of the form 0(r) = exp(— |r^) J2i CiH4i+i{r). 
In such a basis, each term is form invariant with respect to transformations 
to momentum space. The Hamiltonian H has earlier been studied by Boukraa 
and Basdevant [10] with the aid of special methods for solving problems in 
momentum space. Thus we know that the bottom of the spectrum of H to 
4 places is = 2.2322. By use of the tangential bound we obtain an upper 
family of operators of the form H — sp^ -\- l/(4s) -|- r. The corresponding low- 
est eigenvalue is given by S{s) = ss^—zq) + l/(4s), where Ai(2;o) = 0, and 
Zq ~ —2.3381074. By minimizing over s we find the best upper bound to be 

E < (4/3) (3|^o|V4)^ ~ 2.3461, that is to say about 5.1% high. 

We shall return to these examples in Sec. 4 and find sharper upper estimates. 
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3 Difference Schrodinger upper bounds 



The upper bound we shall discuss was discovered in connection with our stud- 
ies of the semirelativistic many-body problem. For the 1-particlc case the 
bound is most easily constructed by means of the following defining equa- 
tions: 



H = ^w? +p'^ + V{r) = Hi- H2, (3.1) 

where 



Hi = \J m? + + ap^ + br^ = Hi = bp^ + \Jvn? -|- -|- ar^, 

H2 = ap" ^hr^ -V{r), (3.2) 

and the parameters a and b are positive. We shall assume that the harmonic 
oscillator potential br^ dominates the potential \^(r) for large r. In this case the 
operators Hi and H2 are both Schrodinger operators whose spectral bottoms 
we write respectively as -E'i(a, 6) and £^2(0, b). These energies can of course be 
found equivalently from eigenproblems expressed in coordinate or momentum 
space. We let E be the bottom of the spectrum of H and we express the 
relation between the Hamiltonian operators in the form 

Hi = H + H2. (3.3) 

It is clear by an elementary variational argument applied to (3.3) that we can 
conclude the Weyl energy inequality [7]- [9] 

Ei{a,b)>E + E2ia,b). (3.4) 

By re-writing (3.4) and optimizing with respect to the free positive parameters 
a and b, we find that our best such difference upper bound to E is given by 

E <Eu^ min [Ei{a, b) - E2{a, b)] . (3.5a) 

{a, b} 

By adding and subtracting the oscillator ap'^ + br^ in the reverse way we arrive, 
by exactly similar reasoning, at an alternative difference upper-bound formula 
given by 

E<Ei-^ = min b) - E[-\a, 6)1 , {3.5b) 

{a, b} L J 

where the corresponding operators H^ ^ and H2 ^ are defined by 



Hi = —y TTT? + p^ + ap^ + br^ = H^ — bp^ — ^Jvn? -|- -|- 



ar 



H^'^ ^ap^ + br^ + V{r). (3.6) 
Equations (3.5a) and (3.5b) summarize the principal results of this paper. 
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4 Examples 



We now consider the three problems mentioned in Sec. 2. In each case we must 
solve the corresponding Schrodinger problems, Hi and H2, defined in (3.2), 
and then minimize the corresponding eigenvalue difference Ei{a, b) — £'2(0, b) 
with respect to the parameters a and b. 

4-1 The ultmrelativistic harmonic oscillator H — p-\-r^ (m = 0) 

The corresponding pair of Schrodinger operators given by (3.2) become 

^1 = 6p2 _^ ar^ + r, H2 = ap'^ + {b - l)r'^ . (4.1) 
We find from (3.5a) 

2.3381074 f« E < 2.3433 = 5.63456 - 3.29126 (a = 0.59, 6 = 3.04). (4.2) 

4-2 The semirelativistic harmonic oscillator H — y/1 + (m = 1) 

The corresponding pair of Schrodinger operators are 

Hi = bp^ + ar^ + Vl + r^, H2 = ap^ + (6 - l)r^ (4.3) 
Prom (3.5a) we find 

2.6640167 fti E < 2.6689 = 6.33418 - 3.66528 (a = 0.59, 6 = 3.53). (4.4) 

4-3 The ultrarelativistic linear potential H — p + r {m — 0) 

The corresponding pair of Schrodinger operators are given by 

Hi = bp + ar + r, H2 = ap + or — r. (4.5) 

We can derive the best upper bound provided by Eq. (3.5) analytically in this 
case. We find 

E<E^= lim [Ei{a, a) - ^2(0, «)] = ^^^7^ = ^ ~ 2.25676, (4.6) 

(0, 0) VTT 

where 0(r) = exp(— |r^). We can see this by the following argument. If we let 
the bottom of the spectrum of the perturbed oscillator + + Ar be e(A), 
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and we write a — and b — t"^, then by scaling arguments we obtain the 
equation 



Eiia, h) - E2ia, b) = sH^ [e [l/{sH)) - e [-l/{th))\ . (4.7) 

This expression provides an upper bound for every choice of s and t. The 
difference will be small when both the expressions for A are small, that is to 
say, when s and t are large. In the limit of small A, we know by perturbation 
theory that the approximation e(A) ~ 3 + (2/-y/7r)A is asymptotically exact. 
Thus we find in this small- A limit that 

E,{a, h) - E^ia, b) ^ ^ (- + ^) > ^. (4.8) 

The minimum implies that s — t, and the small- A limit implies that s — > oo. 
Thus the best upper bound provided by the smallest spectral difference is 
given by the right-hand side of (4.8), as claimed above. 

It is evident that the difference upper bound leads to more accurate results 
for these problems than does the envelope upper bound. We note that the 
bounds provided by the alternative difference inequality (3.5b) are very similar 
in numerical quality. 



5 Conclusion 



The main attraction of the Salpeter Hamiltonian H = \/?7T?~+~p^ + V{r) is 
that it captures some relativistic features whilst remaining a relatively simple 
operator. By simple we mean that for many potentials, its spectrum can be 
defined variationally. Thus it is in principle straightforward to find energy up- 
per bounds by exploring a finite- dimensional trial space. The main technical 
difficulty concerning the Hamiltonian is that, apart from the harmonic oscil- 
lator V{r) = r^, the Hamiltonian is in general non-local. In the present paper 
we explore a new class of Schrodinger operator differences that provide upper 
bounds. The ultrarelativistic linear problem H = p + r shows that in some 
cases we may expect to obtain analytical results from these bounds. 
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